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Abstract 

Adiabatic nucleon-nucleon potentials are calculated in a six-quark nonrel- 
ativistic chiral constituent quark model where the Hamiltonian contains a 
linear confinement and a pseudoscalar meson (Goldstone boson) exchange in- 
teraction between quarks. Calculations are performed both in a cluster model 
and a molecular orbital basis, through coupled channels. In both cases the po- 
tentials present an important hard core at short distances, explained through 
the dominance of the [51] fs configuration, but do not exhibit an attractive 
pocket. We add a scalar meson exchange interaction and show how it can 
account for some middle-range attraction. 

I. INTRODUCTION 

There have been many attempts to study the nucleon-nucleon interaction starting from 
a system of six interacting quarks described by a constituent quark model. These models 
explain the short range repulsion as due to the colour magnetic part of the one gluon 
exchange (OGE) interaction between quarks and due to quark interchanges between two 3q 
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clusters ]T],|2[]. To the OGE interaction it was necessary to add a scalar and a pseudoscalar 
meson exchange interaction between quarks of different 3q clusters in order to explain the 
intermediate- and long-range attraction between two nucleons |3|-|5| . 

In a previous work || we have calculated the nucleon-nucleon (NN) interaction potential 
at zero-separation distance between two three-quark clusters in the frame of a constituent 
quark model |7|-|9| where the quarks interact via pseudoscalar meson i.e. Goldstone boson 
exchange (GBE) instead of OGE. An important motivation in using the GBE model is 
that it describes well the baryon spectra. In particular, it correctly reproduces the order of 
positive and negative parity states both for nonstrange || and strange |§ baryons where 
the OGE model has failed. 

The underlying symmetry of the GBE model is related to the flavour-spin SUp(3) x 
SUs{2) group. Combining it with the S3 symmetry, a thorough analysis performed for the 
L = 1 baryons has shown that the chiral quark picture leads to more satisfactory fits to 



the observed baryon spectrum than the OGE models. 

The one-pion exchange potential between quarks appears naturally as an iteration of 
the instanton induced interaction in the t-channel The meson exchange picture is also 
supported by explicit QCD latice calculations JT2I . 



Another motivation in using the GBE model is that the exchange interaction contains 
the basic ingredients required by the AW problem. Its long-range part, required to provide 
the long-range iViV interaction, is a Yukawa-type potential depending on the mass of the 
exchange meson, Its short-range part, of opposite sign to the long-range one, is mainly 
responsible for the good description of the baryon spectra |7|-|9| and also induces a short- 
range repulsion in the iViV system, both in the 3 Si and 1 S'o channels [[13[]. The present 
study is an extention of and we calculate here the interaction potential between two 
3q clusters as a function of Z, the separation distance between the centres of the clusters. 
This separation distance is a good approximation of the Jacobi relative coordinate between 
the two clusters. Under this assumption, here we calculate the interaction potential in the 
adiabatic (Born-Oppenheimer) approximation, as explained below. 



A common issue in solving the NN problem is the construction of adequate six-quark 
basis states. The usual choice is a cluster model basis |I],0,|14]] . In calculating the potential at 
zero-separation distance, in Ref. || we used molecular-type orbitals [IS] and compared the 
results with those based on cluster model single-particle states. The molecular orbitals have 
the proper axially and reflectionally symmetries and can be constructed from appropriate 
combinations of two-centre Gaussians. At zero-separation between the 3q clusters the six- 
quark states obtained from such orbitals contain certain p n s 6 ~ n configurations which are 
missing in the cluster model basis. By using molecular orbitals, in Ref. || we found that 
the height of the repulsion reduces by about 22 % and 25 % in the 3 S\ and 1 5'o channels 
respectively with respect to cluster model results. It is therefore useful to analyse the role 
of molecular orbitals at distances Z ^ 0. By construction, at Z — > oo the molecular orbital 
states are simple parity conserving linear combinations of cluster model states. Their role is 
expected to be important at short range at least. They also have the advantage of forming an 
orthogonal and complete basis while the cluster model (two-centre) states are not orthogonal 
and are overcomplete. For this reason we found that in practice they are more convenient 
to be used than the cluster model basis, where one must carefully [frjj consider the limit 
Z — > 0. Here too, for the purpose of comparison we perform calculations both in the cluster 
model and the molecular orbital basis. 

In Sec. 2 we recall the procedure of constructing molecular orbital single-particle states 
starting from the two-centre Gaussians used in the cluster model calculations. In Sec. 3 
the GBE Hamiltonian is presented. Sec. 4 is devoted to the results obtained for the NN 
potential. In Sec. 5 we introduce a middle range attraction through a scalar meson exchange 
interaction between quarks parametrized consistently with the pseudoscalar meson exchange. 
The last section is devoted to a summary and conclusions. 
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II. THE SINGLE-PARTICLE ORBITALS 



In the cluster model one can define states which in the limit of large intercluster separa- 
tion Z are right R and left L states 

R = ^^r-^j and L = ip ff + ^ j . (1) 

In the simplest cluster model basis these are ground state harmonic oscillator wave functions 
centered at Z/2 and —Z/2 respectively. They contain a parameter f3 which is fixed varia- 
tionally to minimize the nucleon mass described as a 3q cluster within a given Hamiltonian. 
The states (1) are normalized but are not orthogonal at finite Z. They have good parity 
about their centers but not about their common center f = 0. 

From R and L one constructs six-quark states of given orbital symmetry [f]o- The 
totally antisymmetric six-quark states also contain a flavour-spin part of symmetry [f]FS 
and a colour part of symmetry [222]c. In the cluster model the most important basis states 
|T3| for the Hamiltonian described in the following section are 



R 3 L 3 [6] G [33]^) (2) 
R 3 L 3 [42] G [33U) (3) 
R 3 L 3 [42] G \51) Fa ) (4) 



R 3 L 3 [42] [411] FS ) (5) 



Harvey |TJj] has shown that with a proper normalization the symmetry [6]o contains only s 6 
and [42] o only s 4 p 2 configurations in the limit Z — > 0. 



According to Ref. [T^] let us consider now molecular orbital single-particle states. Most 
generally these are eigenstates of a Hamiltonian H having axial and reflectional symmetries 
characteristic to the NN problem. These eigenstates have therefore good parity and good 
angular momentum projection. As in the cluster model basis where one uses the two lowest 



states R and L, in the molecular orbital basis we also consider the two lowest states, a, of 
positive parity and tt, of negative parity. From these we can construct pseudo-right r and 
pseudo-left I states as 



2~ 1/2 (a ± tt) for all Z, 



(6) 



where 



< r\r >=< l\l >= 1, < r\l >= 0. 



(7) 



In principle one can obtain molecular orbital single particle states from mean field calcula- 
tions (see for example [0). Here we approximate them by good parity, orthonormal states 
constructed from the cluster model states (1) as 



[2(1± < R\L >)}- 1 ' 2 {R±L) 1 



(8) 



Such molecular orbitals are a very good approximation to the exact eigenstates of a "two- 
centre" oscillator frequently used in nuclear physics or occasionally [[T7[| in the calculation of 
the NN potential. They provide a convenient basis for the first step calculations based on 
the adiabatic approximation as described below. 
Introduced in (6) they give 



R + L 



± 



R-L 



1+ < R\L >)Va (1- < R\L >)Va 



(9) 



At Z —>■ one has a — > s and tt —>■ p (with m = 0, ±1) where s and p are harmonic oscillator 
states. Thus in the limit Z — > one has 



2 1/2 (s±p), 



(10) 



and at Z — > oo one recovers the cluster model basis because r — > i? and £ —>■ L. 

Equation (9) with i? and L defined by (1) ensures that the same Z is used both in the 
molecular and the cluster model basis. 



From (r, I) as well as from (cr, n) orbitals one can construct six-quark states of required 
permutation symmetry. For the Se symmetries relevant for the NN problem the transfor- 
mations between six-quark states expressed in terms of (r, I) and (a, ti) states are given in 



Table I of Ref . [To] . This table shows that in the limit Z — ► six-quark states obtained from 
molecular orbitals contain configurations of type s n p 6 ~ n with n = 0,1,. ..,6. For example 
the [6]o state contains s 6 , s 6 p 4 , s 2 p 4 and p 6 configurations and the [42] o state associated 
to the S"-channel contains s 4 p 2 and s 2 p 4 configurations. This is in contrast to the cluster 
model basis where [6]o contains only s 6 and [42] o only s p 2 configurations, as mentioned 
above. This suggests that the six-quark basis states constructed from molecular orbitals 
form a richer basis without introducing more single particle states. 



Using Table I of Ref. [15] we find that the six-quark basis states needed for the Si or 



1 So channels are: 

|33[6] [33W 
|33[42] [33] FS ) 

|33[42] [51W 
|33[42] G [411] F5 ) 

42+[6] [33] FS 

42+[42] [33] FS 

42+[42] [51] FS 

42+[42] [411] FS 
51+[6] [33] FS 



1 

4 



/I 

/I 

1 



VE (a 6 - vr 6 ) - Vz ( 

[42] [33] F .,. 



^4 2 ^2 4 

cr 7r — a n 



[6] G [33] 



FS / ' 



^4 2 



^2 4 

cr ti 



~-4 2 

a ti 



^-2 4 



^4 2 2 4 

(J 7T — (J 7T 



-2^4 
7T 



[6] f 



o\y"\FS / ' 



[42]o[5lU/- 
[42] [411] F5 ), 

15 (a 6 + vr 6 ) - (a 4 vr 2 + a 

[42] [33] F5 ), 

[42] [51U), 

[42] [411] F5 ), 
(a 6 - vr 6 ) + (aV - aV)] [6] G [33] 



^4 2 , ^.2 4 
O 71 + (7 7T 



^4 2 , ^2 4 
<J 71 + <T 71 



^-4 2 
CT 7T 



^-2 4 



FS / ' 



(11) 

(12) 
(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 



where the notation 33 and mn + in the left-hand side of each equality above means r 3 £ 3 and 



r m £ n + r n £ m as in Ref. [15]. Each wave function contains an orbital part (O) and a flavour- 



spin part (-FS 1 ) which combined with the colour singlet [222] q state gives rise to a totally 



antisymmetric state. We restricted the flavour-spin states to [33] fs; [51]fs and [411]fs as 
for the cluster model basis (2-5). 

As explained above, besides being poorer in s n p 6 ~ n configurations, the number of basis 
states is smaller in the cluster model although we deal with the same [f]o and [/]f5 sym- 
metries and the same harmonic oscillator states s and p in both cases. This is due to the 
existence of three-quark clusters only in the cluster model states, while the molecular basis 
also allows configurations with five quarks to the left and one to the right, or vice versa, or 
four quarks to the left and two to the right or vice versa. At large separations these states 
act as "hidden colour" states but at short- and medium-range separation distances they are 
expected to bring a significant contribution, as we shall see below. The "hidden colour" are 
states where a 3q cluster in an s 3 configuration is a colour octet, in contrast to the nucleon 
which is a colour singlet. Their role is important at short separations but it vanishes at 



large ones (see e.g. fL"4"|). 



III. HAMILTONIAN 

The GBE Hamiltonian considered in this study has the form H|J : 

# = E^ + E|^- + E VaUm) + E ^to) , (20) 

i i * * %<j i<j 

with the linear confining interaction : 

Konf(%) = ~X' • \< (V + C rij ), (21) 

and the spin-spin component of the GBE interaction in its SUp(3) form : 

If=i 

+ £ V K (r tJ )XfXf + ^(r^AfAj + VV(r«)A?A5) ff t ■ 9 jy (22) 



with A = y 2/3 1, where 1 is the 3x3 unit matrix. The interaction (22) contains 7 = tt, K, 77 
and rf meson-exchange terms and the form of V 1 (r^) is given as the sum of two distinct 
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contributions : a Yukawa-type potential containing the mass of the exchanged meson and a 
short-range contribution of opposite sign, the role of which is crucial in baryon spectroscopy. 
In the parametrization of Ref. H the exchange potential due to a meson 7 has the form 

a 2 1 P -Vi r A 
V i( r ) = 7T7^r^M r - r o)/4— ^« 3 exp(-a 2 (r-r ) 2 )}. (23) 



4-7T Ylmirrij 7 r y^r 

The shifted Gaussian of Eq. (23) results from a pure phenomenological fit (see below) of 
the baryon spectrum with 

r = 0.43 fm, a = 2.91 frrT 1 , (24) 

For a system of u and d quarks only, as it is the case here, the i^-exchange does not 
contribute. The apriori determined parameters of the GBE model are the masses 

m u , d = 340 MeV, fi n = 139 MeV , ^ = 547 MeV , /v = 958 MeV - (25) 

The other parameters are given in Table I. 

It is useful to comment on Eq. (23). The coupling of pseudoscalar mesons to quarks 
(or nucleons) gives rise to a two-body interaction potential which contains a Yukawa-type 
term and a contact term of opposite sign (see e.g. The second term of (23) stems from 



the contact term, regularized with parameters fixed phenomenologically. Certainly more 
fundamental studies are required to understand this second term and attempts are being 
made in this direction. The instanton liquid model of the vacuum (for a review see |L9| ) 
implies point-like quark-quark interactions. To obtain a realistic description of the hyperfine 
interaction this interaction has to be iterated in the t-channel The t-channel iteration 



admits a meson exchange interpretation p0 



In principle it would be better to use a parametrization of the GBE interaction as given 
in [2]J based on a semirelativistic Hamiltonian. However, in applying the quark cluster 
approach to two-baryon systems we are restricted to use a nonrelativistic kinematics and an 
s 3 wave function for the ground state baryon. 

The matrix elements of the Hamiltonian (20) are calculated in the bases (2)-(5) and (11)- 
(19) by using the fractional parentage technique described in Refs. Hl4| , p2| and also applied 



S 



in Ref. ||13|| . A programme based on Mathematica |23j has been created for this purpose. In 
this way every six-body matrix element reduces to a linear combination of two-body matrix 
elements of either symmetric or antisymmetric states for which Eqs. (3.3) of Ref. can be 
used to integrate in the flavour-spin space. 

IV. RESULTS 

We diagonalize the Hamiltonian (20)- (25) in the six-quark cluster model basis (2)- (5) 
and in the six-quark molecular orbital basis (11)-(19) for values of the separation distance 
Z up to 2.5 fm. Using in each case the lowest eigenvalue, denoted by (H) z we define the 
NN interaction potential in the adiabatic (Born-Oppenheimer) approximation as 

V NN (Z) = (H) z - 2m N - K rel (26) 

Here uin is the nucleon mass obtained as a variational s 3 solution for a 2>q system described by 



where 



the Hamiltonian (20). The wavefunction has the form oc exp — (p 2 + A 2 ) /2[3 2 
p = (r*i — r 2 ) / \/2 and A = (fi + r 2 — 2f 3 ) / V6. The variational solution for = (H) Sq and 
the corresponding (3 is given in Table II. The same value of /3 is also used for the 6q system. 
This is equivalent with imposing the "stability condition" which is of crucial importance in 
resonating group method (RGM) calculations |0]|| . The quantity K rd represents the relative 
kinetic energy of two 3q clusters separated at infinity 

3fr 2 

Krd = 4^ (2?) 

where m above and in the following designates the mass of the u or d quark. For the value 
of (3 of Table II this gives K Tel = 0.448 GeV. 

A. Cluster model 

In Fig. 1 we present the expectation value of the kinetic energy (KE) as a function of Z. 



One can see that for the state 



R 3 L 3 [42] Q ) it decreases with Z but for the state i2 3 L 3 [6] D 



?3r3r 



it first reaches a minimum at around Z = 0.85 fm and then tends to an asymptotic value 
equal to its value at the origin due to its s 6 structure. This value is 

(KE) Z=0 = (KE) z=OQ = ^huj (28) 

where hcu = h 2 /m(3 2 . Actually this is also the asymptotic value for all states. 

The diagonal matrix elements of the confinement potential are presented in Fig. 2. 
Beyond Z > 1.5 fm one can notice a linear increase except for the R 3 L 3 [42] Q [5 1]^) state 
where it reaches a plateau of 0.3905 GeV. 

As an example the diagonal matrix elements of the chiral interaction V x are exhibited 
in Fig. 3 for S — 1, 7 = 0. At Z = one recovers the values obtained in Ref. ||. 
At Z — > oo the symmetries corresponding to baryon-baryon channels, namely [51]^ and 
[33] fs, must appear with proper coefficients, as given by Eq. (29). The contribution due 
to these symmetries must be identical to the contribution of V x to two nucleon masses also 
calculated with the Hamiltonian (20). This is indeed the case. In the total Hamiltonian 
the contribution of the [411] ^3 V x state tends to infinity when Z — > 00. Then this state 
decouples from the rest which is natural because it does not correspond to an asymptotic 
baryon-baryon channel. It plays a role at small Z but at large Z its amplitude in the NN 
wavefunction vanishes, similarly to the "hidden colour" states. Actually, in diagonalizing 
the total Hamiltonian in the basis (2)- (5) we obtain an NN wavefunction which in the limit 



Z — > 00 becomes [14| 



i>NN = I |[6] G [33] FS ) + I I [42] [33] F5 > - \ | [42] [51] FS > (29) 

The adiabatic potential drawn in Figs. 6 and 7 is defined according to Eq. (26) where (H)z 
is the lowest eigenvalue resulting from the diagonalization. Fig. 6 corresponds to S = 1, / 
= and Fig. 7 to S = 0, 7=1. Note that from these curves one should subtract K re i of Eq. 
(27) in order to obtain the asymptotic value zero for the potential. One can see that the 
potential is repulsive at any Z in both sectors. 

Our cluster model results can be compared to previous literature based on OGE models. 
A typical example for the S S± and 1 S'o adiabatic potentials can be found in Ref. |24 ]. The 
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results are similar to ours. There is a repulsive core but no attractive pocket. However, in 
our case, in either bases, the core is about twice higher at Z = and about 0.5 fm wider 
than in ||24l . 



B. Molecular orbital basis 



In the molecular basis the diagonal matrix elements of the kinetic energy are sim- 
ilar to each other as decreasing functions of Z. As an illustration in Fig. 4 we show 
(K.E.) corresponding to 1 33[6] G [33] FS ) and to the most dominant state at Z — 0, namely 
42 + [42] [51] F5 ^> (see ||). The kinetic energy of the latter is larger than that of the former 
because of the presence of the configuration s 2 p 4 with 50 % probability while in the first 
state this probability is smaller as well as that of the p 6 configuration, see eqs. (11) and 
(17). The large kinetic energy of the state (17) is compensated by large negative values of 
(V x ) so that this state becomes dominant at small Z in agreement with Ref. ||. 

The expectation values of the confinement potential increase with Z becoming linear 
beyond Z > 1.5 fm except for the state 1 33[42] [51] FS ) which gives a result very much 



can be understood through the details given in the Appendix. Due to the similarity to the 
cluster model results we do not show here (V con f) explicitly for the molecular orbital basis. 
The expectation value of the chiral interaction either decreases or increases with Z de- 



fer S = 1, I = and S = 0, I = 1 sectors. This state is the dominant component of tpNN at 
Z = with a probability of 87 % for SI = (10) and 93 % for SI = (01) f§. With increasing 
Z these probabilities decrease and tend to zero at Z — > oo. In fact in the molecular orbital 
basis the asymptotic form of ipNN is also given by Eq. (29) inasmuch as r — > R and I — > L 
as indicated below Eq. (10). 

Adding together these contributions we diagonalize the Hamiltonian and use its lowest 
eigenvalue to obtain the A^A^ potential according to the definition (26). The S = 1, / = 



similar to the cluster model state 




drawn in Fig. 2. Such a behaviour 



pending on the state. In Fig. 5 we illustrate the case of the 42 + [42] [51] F5 \ state, both 
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and S = 0, I = 1 cases are illustrated in Figs. 6 and 7 respectively, for a comparison with 
the cluster model basis. As shown in Ref. at Z = the repulsion reduces by about 22 
% and 25 % in the 3 Si and 1 5'o channels respectively when passing from the cluster model 
basis to the molecular orbital basis. From Figs. 6 and 7 one can see that the molecular 
orbital basis has an important effect up to about Z pa 1.5 fm giving a lower potential at 
small values of Z. For Z pa 1 fm it gives a potential larger by few tens of MeV than the 
cluster model potential. However there is no attraction at all in either case. 



Actually, by construction, the molecular orbital basis is richer at Z = |15[ than the 
cluster model basis. For this reason, at small Z it leads to a lower potential than the cluster 
model basis. Within a truncated space this property may not hold beyond some value of Z. 
However by an increase of the Hilbert space one can possibly bring the molecular potential 
lower again. In fact we chose the most important configurations from symmetry arguments 
1 13fl based on Casimir operator eigenvalues. These arguments hold if the interaction is the 
same for all quarks in the coordinate space. This is certainly a better approximation for 
Z = than for larger values of Z. So it means that other configurations, which have been 
neglected, may play a role at Z > 0.4 fm. Then, if added, they could possibly lower the 
molecular basis result. 

As defined in Sec. 2 the quantity Z is the separation distance between two 3q clusters. It 
represents the Jacobi relative coordinate between the two nucleons only for large Z. There we 
view it as a generator coordinate and the potential we obtain represents the diagonal kernel 
appearing in the resonating group or the generator coordinate method. The comparison 
given above should then be considered in the context of the generator coordinate method 
wich will be developped in further studies and will lead to nonlocal potentials. However the 
adiabatic potentials, here obtained in the two bases can be compared with each other in an 
independent and different way. On can introduce the quadrupole moment of the six-quark 
system 

g 2 o = Er J 2 F 20 (r J ) (30) 



i=i 
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and treat the square root of its expectation value 

(Q) = (iPnn^oI^nn) (31) 
as a collective coordinate describing the separation between the two nucleons. Obviously 



(Q) -> Z for large Z. 



In Fig. 8 we plot J (Q) as a function of Z. The results are practically identical for IS 



(01) and IS = (10). Note that y/(Q) is normalized such as to be identical to Z at large Z. 



One can see that the cluster model gives y (Q) = at Z = 0, consistent with the spherical 



symmetry of the system, while the molecular basis result is y (Q) = 0.573 fm at Z — 0, 
wich suggests that the system acquires a small deformation in the molecular basis. This also 
means that its r.m.s. radius is larger in the molecular basis. 



In Figs. 9 and 10 we plot the adiabatic potentials as a function of y (Q) instead of Z, for 



IS = (01) and (10) respectively. As y/(Q) 7^ at any Z in the molecular orbital basis, the 
corresponding potential is shifted to the right and appears above the cluster model potential 



at finite values of y (Q) but tends asymptotically to the same value. The comparison made 
in Figs. 9 and 10 is meaningful in the context of a Schrodinger type equation where the local 
adiabatic potential appears in conjunction with an "effective mass" depending on \J (Q) also. 
However an effective mass can be obtain through RGM calculations and our future plan is 
to perform such calculations. 



V. THE MIDDLE RANGE ATTRACTION 

In principle we expected some attraction at large Z due to the presence of the Yukawa 
potential tail in Eq. (23). To see the net contribution of this part of the quark-quark inter- 
action we repeated the calculations in the molecular orbital basis by completely removing 
the first term - the Yukawa potential part - in Eq. (23). The result is shown in Fig. 11. for 
SI— (10). One can see that beyond Z 1.3 fm the contribution of the Yukawa potential 
tail is very small, of the order of 1-2 MeV. At small values of Z the Yukawa part of (23) 
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contributes to increase the adiabatic potential because it diminishes the attraction in the 
two body matrix elements. 

The missing medium- and long-range attraction can in principle be simulated in a simple 
phenomenological way. For example, in Ref. |]J this has been achieved at the baryon level. 
Here we adopt a more consistent procedure assuming that besides the pseudoscalar meson 
exchange interaction of Sec. Ill there exists an additional scalar, cr-meson exchange inter- 
action between quarks. This is in the spirit of the spontaneous chiral symmetry breaking 
mechanism on which the GBE model is based. The cr-meson is the chiral partner of the pion 
and it should be considered explicitly. 

Actually once the one-pion exchange interaction between quarks is admitted, one can 
inquire about the role of at least two-pion exchanges. Recently it was found 0] that the 
two-pion exchange also plays a significant role in the quark-quark interaction. It enhances 
the effect of the isospin dependent spin-spin component of the one-pion exchange interaction 
and cancels out its tensor component. Apart from that it gives rise to a spin independent 
central component, which averaged over the isospin wave function of the nucleon it produces 
an attractive spin independent interaction. These findings also support the introduction of 
a scalar (cr-meson) exchange interaction between quarks as an approximate description of 
the two-pion exchange loops. 

For consistency with the parametrization || we consider here a scalar quark-quark in- 
teraction of the form 

V a (r) = f^—{0(r - r'M *— - -^a' 3 exp(-c/ 2 (r - r(,) 2 )}. (32) 

where = 675 MeV and r' , a' and the coupling constant g^/kn are arbitrary parameters. In 
order to be effective at medium-range separation between nucleons we expect this interaction 
to have r' ^ r and a' ^ a. Note that the factor 1/mirrij has only been introduced for 
dimensional reasons. 

We first looked at the baryon spectrum with the same variational parameters as before. 
The only modification is a shift of the whole spectrum which would correspond to taking 
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V Q « -60 MeV in Eq. (21). 

For the 6q system we performed calculations in the molecular basis, which is more ap- 
propriate than the cluster model basis. We found that the resulting adiabatic potential is 
practically insensitive to changes in \x a and r' but very sensitive to a'. In Fig. 12 we show 
results for 

r' = 0.86 fm, a' = 1.47 fm~\ g 2 a /4n = g 2 8 /4n (33) 

One can see that V a produces indeed an attractive pocket, deeper for SI = (10) than for 
(01), as it should be for the NN problem. The depth of the attraction depends essentially 
on a' . The precise values of the parameters entering Eq. (32) should be determined in 
further RGM calculations. As mentioned above the Born-Oppenheimer potential is in fact 
the diagonal RGM kernel. It is interesting that an attractive pocket is seen in this kernel 
when a cx-meson exchange interaction is combined with pseudoscalar meson exchange and 
OGE interactions (hybrid model), the whole being fitted to the NN problem |25| . 



VI. SUMMARY 

We have calculated the NN potential in the adiabatic approximation as a function of Z, 
the separation distance between the centres of the two 3q clusters. We used a constituent 
quark model where quarks interact via pseudoscalar meson exchange. The orbital part of 
the six-quark states was constructed either from cluster model or molecular orbital single 
particle states. The latter are more realistic, having the proper axially and reflectionally 
symmetries. Also technically they are more convenient. We explicitly showed that they are 
important at small values of Z. In particular we found that the NN potential obtained in 
the molecular orbital basis has a less repulsive core than the one obtained in the cluster 
model basis. However none of the bases leads to an attractive pocket. We have simulated 
this attraction by introducing a a-meson exchange interaction between quarks. 

To have a better understanding of the two bases we have also calculated the quadrupole 
mement of the 6q system as a function of Z. The results show that in the molecular orbital 
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basis the system acquires some small deformation at Z = 0. As a function of the quadrupole 
moment the adiabatic potential looks more repulsive in the molecular orbital basis then in 
the cluster model basis. In this light one might naively expect that the molecular basis will 
lead to scattering phase-shifts having a more repulsive behaviour. 

The present calculations give us an idea about the size and shape of the hard core 
produced by the GBE interaction. Except for small values of Z the two bases give rather 
similar potentials. Taking Z as a generator coordinate the following step is to perform a 
dynamical study based on the resonating group method which will provide phase-shifts to 
be compared to the experiment. The present results constitute an intermediate step towards 
such a study. 



VII. APPENDIX 

In this appendix we study the behaviour of the confinement potential in the molecular 
orbital basis at large separation distance Z between the centres of two 3q clusters. As an 
example we consider the state |42 + [42] [33]^,s). Through the fractional parentage technique 
|14] , |22[| the six-body matrix elements can be reduced to the calculation of two-body matrix 
elements. Using this technique and integrating in the color space one obtains 

(42 + [42] [33] F5 |V; on/ |42 + [42]o[33] F 5) = ^ [22(mr\V\nn) 
+ 76((T7r|V|(T7r) + 26(<X7r| V\na} 

- 58(nn\V\aa) + 22(aa\V\aa)] (34) 

where the right-hand side contains two-body orbital matrix elements. According to Eq. (8) 
for Z — > oo one has 

\a) - |it! + L), |tt) - i= \R - L) (35) 

Replacing these asymptotic forms in the above equation one obtains matrix elements con- 
taining the states \R) and \L). Most of these martix elements vanish asymptotically. The 
only surviving ones are 
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(R R\V\R R)^a , (R L\V\R L) -> 6 Z (36) 

where a and 6 are some constants. This brings to the following asymptotic behaviour of the 
matrix elements in the right-hand side of (34) 



(aa \V \aa) — > 


(a 


+ 


b Z)/2 


{nn \V 7T7r) — > 


(a 


+ 


b Z)/2 


(<77r| V|<77r) — > 


(a 


+ 


b Z)/2 


(air\V\ira) — > 


(a 




b Z)/2 


(7T7T V ffff) — > 


(a 




b Z)/2 



(37) 



from which it follows that 

(42 + [42]o[33] FS |K „ / |42 + [42]o[33] F s) -> (H a + 19 b Z)/10 (38) 

i.e. this matrix element grows linearly with Z at large Z. In a similar manner one can show 
that the confinement matrix element of the state 1 33[42] G [51] FS ) the coefficient of the term 
linear in Z cancels out so that in this case one obtains a plateau as in Fig. 2. 

Acknowledgements. We are very grateful to L. Wilets, K. Shimizu and L. Glozman 
for useful comments. 
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TABLES 

TABLE I. Parameters of the Hamiltonian (20-25) 



V (MeV) 


C (fm" 2 ) 


08/47T 


9o 2 /4vr 


Reference 





0.474 


0.67 


1.206 


i 



TABLE II. Variational solution of the Hamiltonian (20)-(25) for the nucleon mass ttim with (3 
as a variational parameter (see text) 



P (fm) 


m N (MeV) 


0.437 


969.6 
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FIG. 1. The cluster model basis. The expectation value of the kinetic energy (K.E.) as a 
function of the separation distance Z between two 3g clusters. The asymptotic value of 2.242 GeV 
, given by Eq. (28) is indicated. The full line corresponds to \ [G]o) and the dashed line to |[42]o) 
states. 
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0.39 



0.3895 



0.389 



0.3885 



0.388 



0.3875 



0.387 



Z(fm) 



FIG. 2. The cluster model basis. The expectation value of V con f of Eq. (21). The correspond- 
ing states are: (1) - |[6]o[33] FS >, (2) - |[42] [33] FS ), (3) - |[42] [51] F5 ), (4) - |[411] [51] F5 )- Note 
that for curve (3) the scale is on the rhs vertical line. 
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FIG. 3. The cluster model basis. The expectation value of the chiral interaction , Eqs. 
(22)-(25), for S = 1, / = 0. The curves are numbered as in Fig. 2 and the scale for (3) is also on 
the rhs vertical line. 
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FIG. 4. The molecular orbital basis. The expectation value of the kinetic energy (K.E.) for 
the |[6] [33]fs) (full curve) and |42+[42] [51]fs) (dashed curve) states (see Eqs. (11) and (17) 
respectively). The latter is the most dominant state at Z = (see text). 
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FIG. 5. The molecular orbital basis. The expectation value of the chiral interaction, eqs. 
(22)-(25), for |42 + [42]o[51]f5) which is the most dominant state at Z = 0. The dashed curve 
corresponds to S = 1, I = and the full curve to S = 0, I = 1. 
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FIG. 7. Same as Fig. 6 but for S = 0, I = 1. 
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FIG. 8. y (Q) as a function of Z with (Q) denned by Eq. (31) and normalized as indicated in 
the text. The full line corresponds to the cluster model basis and the dashed line to the molecular 
orbital basis. 
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FIG. 10. Same as Fig. 9 but for S = 0, / = 1. 
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FIG. 11. The adiabatic potential in the molecular orbital basis for SI = (10). The solid curve 
is the same as in Fig. 6. The dashed curve is the result obtained by removing the Yukawa part of 
the quark-quark interaction (23). 
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FIG. 12. The adiabatic potential in the molecular orbital basis for SI = (10) (full curve) and 
SI = (01) (dashed curve) with pseudoscalar + scalar quark-quark interaction. 
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